(ZfJ$ 

<1  '?  *f2> 


SPACE  SCIENCES  LABORATORY 


AEROPHYSICS  SECTION 


RADIANCE  OF  SPECIES  IN  HIGH  TEMPERATURE  AIR 


/7>\ftr 


(f)\^ 


<S> 

by 


Maria  Nardone,  R.  G.  Breene, 
Saydean  Zeldin*  T.  R.  Riethof  * 


O  -w  UV\  ^“3ry 


Work  performed  under  United  States  Air  Force 


w  or k  periormea  unaer  u 
C)  S^Contract  AF  04I694X-222 

_  A  4=^- 

CE>  **P. 


Consultant,  "Physical  Studies,  Inc.  " 


r\  * 

(  /  V/J  ***"' 

(J  A/*  *  R63SD3 


June,  1963 


HV 


; '  <  i  \  t 

...  ..*-.*"* 


MISSILE  AND  SPACE  DIVISION 


GENERAL®  ELECTRIC 


ABSTRACT 


The  spectral  and  total  radiance  of  the  radiating  systems  in  high 

6* 

temperature  air  have  been  computed  for  temperatures  between  3000  K 
and  25000  K,  and  for  relative  densities  between  Kj  and  10. 

The  total  and  spectral  radiance  of  equilibrium  air  over  the  same  tem¬ 
perature  and  density  range  has  been  computed  from  the  composition 
of  equilibrium  air  and  the  absorption  coefficients  of  the  individual  systems. 
The  results  are  compared  with  those  obtained  by  pervious  investigators. 


v  The  results  are  presented  in  graphical  form  and  equations  approxi¬ 
mating  the  results  in  regions  of  low  emissivity  are  derived. 
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INTRODUCTION 


This  report  considers  the  theories  utilized  in  the  calculations  of  the 
equilibrium  spectral  radiance  of  the  chemical  species  in  high  temperature 
air.  The  knowledge  of  the  radiance  of  each  of  the  radiating  systems  com¬ 
bined  with  a  knowledge  of  the  composition  of  air  as  a  function  of  temperature 
and  density  permits  the  calculation  of  the  total  equilibrium  radiation  of  high 
temperature  air.  These  calculations  apply  equally  well  to  air  which  is  at 
equilibrium,  or  to  air  which  is  in  chemical  non- equilibrium,  provided  only 
that  all  other  degrees  of  freedom  are  assumed  in  equilibrium,  and  provided 
that  the  state  (temperature,  species  concentrations)  of  the  non- equilibrium 
air  is  defined. 


The  input  parameters  involved  in  using  the  theory  in  the  calculation  of 
spectral  absorption  coefficients  are  discussed  in  some  detail.  The  best 
present  values  are  cited,  giving  the  methods  by  which  they  were  obtained 
when  pertinent.  The  influence  of  one  particular  parameter,  the  oscillator 

strength,  is  shown  by  comparison  of  the  results  obtained  by  other  workers 

n  ..  30  .42,43,47,49,50,51,52 

(Kivel  and  Baxley  ,  Treanor,  Wurster,  et  al  , 

22 

Meyerott,  Sokoloff  and  Nicholls  )  with  results  from  air  computations  using 
the  currently  accepted  oscillator  strength  values.  It  is  particularly  significant 

to  note  that  the  curves  for  the  radiance  of  a  given  system  as  obtained  by 

30  7 

Kivel  and  Bailey  and  as  obtained  by  the  later  computations  differ  mostly 

by  a  constant  factor,  showing  that  the  main  change  is  due  to  the  better 


experimental  oscillator  strengths  values  now  available. 


« 
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This  report  includes  a  series  of  figures  which  give: 

1.  The  toted  Radiance  of  Equilibrium  Air  as  a  function  of 
relative  density  for  temperatures  between  3000°K  and 

25,  000°K.  A  cross-plot,  giving  the  radiance  as  a  function 
of  temperature  at  different  relative  densities  is  also  presented. 

2.  The  spectral  Radiance  of  Equilibrium  Air. 

3.  The  total  Radiance  of  Equilibrium  Air  as  a  xunction  of  tem¬ 
perature,  indicating  the  relative  contributions  of  the  various 
radiating  systems. 

4.  The  Radiance  of  each  radiating  system  as  a  function  of  optical 
depth  for  a  range  of  temperatures.  A  second  set  of  figures 
crossplots  the  data,  giving  the  radiance  of  each  system  as  a 
function  of  temperature  at  several  different  concentrations. 

For  each  system,  the  spectral  radiance  is  indicated  at 
several  temperatures. 

Figures  are  included  to  show  the  comparison  of  the  present  results 
with  earlier  data,  some  of  which  used  incorrect  oscillator  strengths 
and  approximate  calculation  methods. 
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THEORETICAL  BASIS 


In  any  quantitative  description  of  the  radiation  from  a  gas  one  may 
choose  to  deal  with  the  linear  gaseous  absorption  coefficient,  the  gaseous 
emissivity,  or  the  actual  intensity  in  units  of  energy  per  unit  of  surface 
area  emitted  by  some  boundary  of  the  gas.  All  three  will  be  discussed, 
so  that  the  table  user  may  choose  the  one  most  suited  to  his  needs. 


Lambert's  law  for  the  absorption  of  radiant  energy  is  given  as: 

_  ,  -KL  _  -Nff L 

I  =  I  e  =  I  e 
o  o 


(1) 


where  I  is  the  intensity  in  the  incident  beam,  I  is  the  intensity  in  the  beam 
after  it  traverses  L  units  of  the  gas,  "K"  is  the  linear  absorption  coef¬ 
ficient  of  the  gas,  "a"  is  the  cross  section  for  this  absorption,  and  "N" 
is  the  number  density  of  absorbing  particles.  Insofar  as  this  equation  is 
concerned,  one  may  consider  it  as  applying  (1)  at  a  particular  radiant  fre¬ 
quency,  (2)  as  a  mean  over  a  band  of  frequencies,  or  (3)  as  a  mean  over 
the  entire  spectrum.  K  and  O’  must  be  individually  considered  for  each  of 
the  types  of  radiation  absorbing  processes,  and,  if  there  are  several  for  the 
frequency  or  frequency  region  of  interest,  they  must  be  summed. 

In  the  present  treatment  it  will  be  assumed  that  local  thermodynamic 

20  38 

equilibrium  is  maintained,  or,  equivalently,  that  Kirchhoff's  law  ’ 
holds.  Under  these  conditions  the  gaseous  emissivity  may  be  written  as: 

c=l-e‘KL  (2) 

Here  again  the  emissivity  may  be  taken  at  a  particular  frequency  or  as 


I 

I 


averaged  over  some  band. 
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Here  o  is  the  usual  Boltzmann  constant  and  e  is  the  average  emissivity. 
Now  the  problem  is  one  of  obtaining  the  cross  section  or  absorption  coef¬ 
ficient  for  whatever  processes  contribute  to  the  radiati  ■  om  high  tem¬ 
perature  air.  First  one  lists  those  processes  which  are  considered  im¬ 
portant  for  such  radiation  in  the  temperature  range  of  interest: 


1.  Ultraviolet-visible  band  spectra 

a.  The  beta  system  of  NO 

b.  The  gamma  system  of  NO 

c.  The  Schumann-Runge  system  of  oxygen 

d.  The  first  positive  system  of  nitrogen 

e.  The  second  positive  system  of  nitrogen 

f.  The  first  negative  system  of  positive  nitrogen 

2.  Infrared  band  spectra 

a.  The  infrared  spectrum  of  NO 
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3.  Bound  free  continua 


a.  The  bound  free  continuum  of  oxygen 

b.  The  deionization  continuum  of  oxygen  ions 

c.  The  deionization  continuum  of  nitrogen  ions 
4.  Free-free  continua 

a.  The  free  free  continuum  of  oxygen 

b.  The  free  free  continuum  of  nitrogen 

c.  The  free  free  continuum  of  N+,  N++,  0+,  and  0++ 

At  higher  temperatures  other  radiation  processes  become  important,  but 

these  will  suffice  here.  In  the  following  subparagraphs  each  of  these  processes 

and  the  methods  of  treating  them  will  be  discussed. 

1.  Ultraviolet -visible  band  spectra 
32  24 

The  Mayer  -  Goody  model  for  an  electronic -vibrational  band 
was  used  in  developing  the  programs  for  ultraviolet -visible  emis- 
sivity  or  absorption  coefficient  calculations.  In  developing  their 
band  model  these  authors  suppose  a  statistical  distribution  of 
spectral  lines  in  a  particular  region  of  the  spectrum.  Analytical 
expressions  for  the  form  of  the  absorption  coefficient  are  then 
developed  from  this  supposition.  We  have  added  the  assumption 
that  the  half  widths  of  the  spectral  lines  are  very  large.  This  as¬ 
sumption  simply  leads  to  an  effectively  continuous  distribution  of 
radiation  within  the  confines  of  the  electronic -vibrational  band.  The 
bands  possess  the  degraded  shape  typical  of  these  bands.  The  fre¬ 
quency  dependent  emissivity  expression  is  now: 
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(6) 


e  =  1  -  exp  [ -dLmJ 


me 


Considering  first  the  summation  in  equation  6,  one  may  suppose 
that  the  band  system  under  consideration  is  degraded  toward  shorter 
wavelengths;  then,  at  any  particular  wavelength,  there  will  be  a  con¬ 
tribution  to  the  radiation  from  all  band  systems  whose  heads  are 
toward  longer  wavelengths  from  that  in  question.  "Q"  is  the  partition 
function  which  is  taken  as  a  product  of  electronic,  vibration  and  ro¬ 
tation  functions.  The  B  are  the  rotational  constants  ;  the  G  ,  are 

v  el,  v 

the  energies  of  the  electronic-vibrational  term  values;  "g"  is  the 

degeneracy  of  the  lower  state,  a  factor  which  is  cancelled  by  the 

electronic  partition  function  when  the  lower  state  is  the  ground 

state;  the  V  ,  are  the  frequencies  of  the  band  heads;  "L"  is  the 
w' 

thickness  of  the  emitting  layer;  N  is  the  number  density  in  particles 
3  23 

per  cm  as  given  by  Browne  »  and  "^nm"  is  the  oscillator  strength 

for  the  band  in  question.  The  oscillator  strengths  are  determined 

37 

from  experimental  measurements  (cf.  Penner  ).  The  uncertainties 
in  the  knowledge  of  the  oscillator  strengths  are  directly  reflected 
in  corresponding  uncertainties  in  the  values  of  the  emissivities. 

2.  Infrared  Band  Systems 

Although  the  physical  processes  involved  are  quite  distinct,  the 


6 


3. 


rotational  constants  really  give  rise  to  the  obvious  differences 
between  an  electronic -vibrational  band  and  a  rotation-vibrational 
band.  In  the  latter  case  the  rotational  constants  are  practically 
the  same  in  the  upper  and  lower  states.  When  one  makes  this  as- 

g 

sumption  one  obtains  a  form  of  Eqs.  6-8  leading  to  the  symmetric 
distribution  of  intensity  in  the  band  typical  of  infrared  spectra.  The 
emissivity  expressions  are: 

e  =  1  -  exp  [  -DE/ii  ]  (9) 


vv 


v  -  V 


4B2  Q 

V 


—  «xP(-Gv/kT)  exp  [-(IM^,)2  4BhCkf]  exp  [-(  V-V^) 

v 


D  =  L 


ne 


N 


[i-e-hWkT] 


he 

Vv1'  2kT 
(10) 

(11) 


me 

The  symbols  are  effectively  the  same  although  no  electronic  transi¬ 
tion  takes  place  here  as  was  the  case  with  the  considerations  of  the 
last  subsection.  Thus,  for  example,  the  oscillator  strength  will 
simply  be  that  of  the  vibrational  transition.  It  should  be  noted  that 
only  one  infrared  spectrum  was  considered,  that  of  NO,  where  the 

g 

oscillator  strengths  used  were  those  computed  by  Breene  and  Todd  , 

based  on  the  integrated  intensity  measurements  of  Weber  and  Penner 

2 

(See  also  Benitez  and  Penner  ) 

Bound  Free  Continua 

a.  Oxygen  free  bound  continuum.  In  the  medium  temperature  re¬ 
gion  (6000°  -  9000°K)  the  important  Bound-Free  continuum  is  that  of 
oxygen.  This  continuum  arises  from  the  attachment  of  an  electron 


7 


to  an  oxygen  atom,  with  the  consequent  emission  of  a  photon.  The 


cross  section  utilized  is  the  one  for  the  inverse  process,  namely 

the  absorption  of  a  photon  by  an  0  ion.  The  cross  section  for  this 

process  can  be  calculated,  but  for  the  present  computations  the 

6 

cross  sections  obtained  by  Branscomb  et  al  in  scattering  experi¬ 
ments,  were  used. 

It  should  be  noted  that  the  cross  section  is  both  temperature 
and  density  independent.  It  is  true  that  the  absorption  coefficient 
will  depend  on  both  these  factors,  but  this  will  only  be  through  the 
dependence  of  the  particle  concentration  on  these  quantities. 

b.  N+  and  Q+  deionization  continua.  The  deionization  cross  - 

+  +  17 

sections  for  N  and  O  have  been  computed  by  Breene  and 

compared  with  the  hydrogenic  cross  sections  for  the  recombination 

process  given  by  Bates  et  al\  In  view  of  the  fact  that  good  agreement 

was  obtained,  the  more  complete  table  of  cross  sections  of  Bates 

et  al  was  applied  to  the  present  problem. 

The  relationship  between  the  cross  section  for  ionization  (<7  ) 

required  in  the  present  calculations,  and  the  recombination  cross 

section  (a  )  given  by  Bates  *  is: 

E  2 
mE  c^ 
k 

a.  =  - 5— ^5 —  a 

A  ,22  e 
h  V 


(12) 


k  - 1 1 

E.  =  [  hcco  +  -T-  (4.  35  x  10  )] 

k  2 


(13) 


60  =  the  frequency  (wavenumbers)  of  the  absorption  edge.  The 

absorption  coefficient  K  is  given  by 

8 


K  =  ctNQ 


(14) 


where  Q  is  the  usual  partition  function,  including  the  appropriate 
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degeneracy  factors.  Details  of  the  procedure  are  given  by  Breene 
4.  The  Free-Free  Continua 

A.  Neutral  Atom-Continua.  Classically,  these  continua  may  be 
considered  as  arising  due  to  the  deflection  of  the  incident  electron 
by  the  atom  with  the  resultant  acceleration  radiation  from  the  de¬ 
flected  electron.  From  the  quantum  point  of  view  such  radiation 
results  from  transitions  on  the  part  of  the  free  electron  between 
different  states  in  the  presence  of  the  particular  atom. 

The  cross  section  for  these  continua  were  calculated,  obtaining 
first  the  wave  functions*^’  ^  for  the  free  electron  in  the  field  of 
the  neutral  atom  in  question.  Originally^’  ^  these  calculations 
were  carried  out  under  the  assumption  of  no  electron  exchange  of 


core  polarization.  The  calculations  for  the  wave  functions  have  b 


een 


17  18 

repeated  *  ,  without  these  limiting  assumptions  and  the  new  re¬ 

sults  are  used  in  the  present  tabulations.  It  should  be  noted  that  the 
17 

new  wave  functions  yield  cross  section  about  one-tenth  those  com- 
12 

puted  earlier  where  free  electron  wave  functions  without  polariza¬ 
tion  were  used.  The  cross  section  values  obtained  for  oxygen  are 

-39  -39  -37 

5.  8  x  10  at  1. 5  microns,  17x10  at  3  microns  and  50  x  10 

at  4.5  microns.  The  results  are  in  excellent  agreement  with  the  data 

r  m  ,  41 

of  Taylor 


9 


B.  Ion  Continua.  The  continua  for  the  ion-electron  deflection 


(Bremsstrahlung)  are  computed  directly  from  the  formula 

45 

originally  developed  by  Kramers,  and  given  by  Unsole  as 


3.  7  x  108  Z2 


where  Z  is  the  charge  on  the  ion, 
2  for  N++  and  0++  etc.). 


(1  for  the  case  of  N 


+ 


and  O 


(15) 

+ 
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INPUT  PARAMETERS 


The  input  parameters  which  were  used  in  the  calculations  are  specified 
in  the  following  in  order  that  the  table  user  may  be  aware  of  the  extent,  if 
any,  to  which  more  accurate  values  for  these  parameters  will  affect  the  tables. 
1.  The  Beta  Bands  of  NO 


This  system  arises  from  the  transition  between  the  two  doublets 
2  2 

X  n3/2  -*  B  n^2  leading  effectively  to  a  double  system.  The 

26 

values  of  the  constants  were  obtained  from  p.  558  of  Herzberg. 


One  can  designate  the  oscillator  strength  for  an  electronic  vibrational 


band  as 


f  =  f  .  f  ,  R  2 
nm  el  vv'  e 


(v) 


(16) 


f  ,  is  proportional  to  the  square  of  the  appropriate  overlap  integral,  this 
being  the  integral  over  the  product  of  the  upper  and  lower  vibrational  state 
wave  functions.  These  integrals  have  been  computed  using  vibrational 
wave  functions  of  various  degrees  of  sophistication,  i.  e. ,  harmonic  oscilla¬ 
tor  functions,  Morse  potential  functions,  and  various  solutions  to  the 
Schroedinger  equation  with  specific  potentials  from  experiment  or  theory. 
For  this  beta  system,  as  well  as  for  systems  to  be  subsequently  described, 
we  have  appealed  to  the  best  available  calculation  of  the  overlap  integrals. 

The  electronic  oscillator  strength  is  not,  of  course,  a  constant  over 
an  electronic  vibrational  band  system.  This  is  perfectly  obvious  from 
the  most  elementary  considerations.  Nevertheless,  it  is  convenient  to 
define  a  constant  quantity  for  a  band  system  which  we  call  the  electronic 
oscillator  strength  and  which  we  refer  to  above  as  f 


I 

I 


11 


This  constant  is 


then  corrected  for  the  various  pairs  of  vibrational  states  leading  to  the 

system  bands.  A  convenient  way  of  doing  this  is  with  the  so-called 

2 

r-centroid  function,  ,  of  Nicholls.  In  this  system,  as  well  as  in 

the  ones  to  be  considered  subsequently,  we  have  applied  the  most  accurate 

r-centroid  available.  In  certain  cases  the  r-centroid  function  has  not 

been  evaluated  from  experiment  for  as  many  vibrational  bands  as  we 

desire  to  treat.  In  these  cases  we  simply  extrapolate  the  r-centroid 

function  as  the  best  available  representation  for  these  bands. 

4 

Bethke  measured  the  absorption  of  nitric  oxide  gas  at  room  tem¬ 
perature  and  calculated  the  oscillator  strength  to  be  0.0015,  a  value  in 

21 

general  agreement  with  the  results  of  Daiber  and  Williams  .  The 

31 

values  of  f  ,  were  computed  by  Kivel,  Mayer  and  Bethe  ,  and  the 

function  R  (r)  was  derived  from  the  experimental  work  and  tabulated  as 
e 

35 

a  function  of  frequency  by  Nicholls 

The  f  .  values  for  2ll„  were  multiplied  by  the  factor  exp  [  -hcp/kT]  V-  120, 
vv1  3/2 

in  order  to  account  for  the  distribution  over  the  lower  members  of  the 
doublet. 

Any  degeneracy  in  the  lower  state  is  cancelled  here  by  the  electronic 
partition  function. 

2.  The  Gamma  Bands  of  NO 

2_+ 

This  system  arises  from  the  transition  between  the  A  E  level  and 
2 

the  X  n3/2  ground  doublet.  The  constants  required  were  obtained 

26 

from  page  558  of  Herzberg. 
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The  value  of  the  electronic  oscillator  strength  was  obtained  by 
4 

Bethke  as  .0024.  A  slightly  higher  value  is  quoted  by  Daiber  and 
21 

Williams.  The  vibrational  overlap  integrals  were  taken  from 

27  2  — 

Jarmain,  Fraser,  and  Nicholls.  The  function  of  R  (r)  was  derived 

e 

39 

by  Robinson  and  Nicholls  . 

The  factor  exp  t-hcp/kT],  V  =  120,  was  again  included.  No  lower 
state  degeneracy  was  considered. 

3.  The  Schumann-Runge  system  of  O 


3  3 

This  system  arises  from  the  transition  X  E  -*  B  £  .  The  constants 

g  g 

2  6 

required  were  obtained  from  Herzberg,  pages  559,  560. 

5 

The  electronic  oscillator  strength  was  obtained  by  Bethke  as  0.  17. 

(Note:  There  is  no  discrepancy  with  the  value  of  0.  048  quoted  by  Treanor 
43 

and  Wurster  ,  since  these  authors  evaluated  their  oscillator  strength 

at  a  different  wave  length).  The  vibrational  overlap  integrals  were 

22  2  — 

obtained  from  Fraser,  Jarmain,  and  Nicholls  ,  while  R  (r)  was  obtained 

e 

34 

from  Nicholls  . 

4.  The  First  Negative  System  of  N  + 


2  2  + 

This  system  arises  from  the  transition  X  S  -*  B  E  .  The 
7  g  M 


constants  were  taken  from  page  554  of  Herzberg. 


26 


The  electronic  oscillator  strength  was  obtained  from  radiative  life-time 

3 

measurements  by  Bennett  and  Dalby  as  0.  04.  The  vibrational  overlap 

27  34 

integrals  were  obtained  from  Jarmain,  Fraser,  and  Nicholls.  Nicholls 


provided  the  function  of  the  r -centroid. 
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5.  The  First  Positive  System  of  N 

b 

This  system  arises  from  the  transition  between  the  two  upper 

3  3  +  26 

states  B  ft  and  A  £  .  The  constants  were  taken  from  Herzberg. 

g  M 

The  electronic  oscillator  strength  of  0.0095  was  obtained  from 

49 

the  measurements  on  high  temperature  nitrogen  of  Wurster.  (The 

49 

oscillator  strength  quoted  by  Wurster  is  .0028.  However,  a  com¬ 
parison  of  the  experimentally  measured  absorption  coefficient  with  our 
calculated  absorption  coefficient  gives  an  electronic  oscillator  strength 
of  .0095.  Through  correspondence  with  Treanor  their  present  electronic 
oscillator  strength  is  .0038.  The  difference  in  these  values  is  due  to 
the  difference  in  the  theoretical  consideration  in  obtaining  the  absorption 
coefficient. ) 

The  electronic  partition  function  involves  the  degeneracy  of  the 

ground  electronic  state,  which  in  this  case  is  unity.  The  lower  state 

involved  in  the  transition,  however,  is  not  the  ground  state  but  an  upper 

state  of  degeneracy  three.  Therefore,  the  degeneracy  is  not  cancelled 

by  the  partition  function  as  in  previous  cases,  and  the  factor  three 

must  be  included.  Here  (g/Qg)  =  3. 

The  vibrational  overlap  integrals  were  obtained  from  Nicholls, 

36 

Jarmain,  and  Fraser  and  the  r-centroid  function  by  Turner  and 
44 

Nicholls, 

25 

Heath  has  noted  that  the  radiation  due  to  the  nitrogen  first 
positive  system  in  air  may  be  less  than  the  values  predicted  from 
measurements  in  pure  nitrogen  since  there  appears  to  be  a  radiation¬ 
less  depopulation o f  the  "B"  state  in  the  presence  of  or  NO. 
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6.  The  Second  Positive  System  of  N_ 

Ld 

This  system  arises  from  the  transition  between  the  two  upper 
3  3 

states  C  II  and  B  II  .  The  parameters  were  obtained  from  page 
^  8 

552  of  Herzberg.^ 

The  electronic  oscillator  strength  was  obtained  from  Bennett 
3 

and  Dalby  as  0.038.  Since  the  lower  state  is  not  the  ground  state, 

the  factor  (g/Q^)  was  taken  as  six. 

The  vibrational  overlap  integrals  were  taken  from  Turner  and 
44 

Nicholls  while  the  r-centroid  function  is  obtained  from  Wallace 
46 

and  Nicholls. 

7.  Comparison  of  Electronic  Oscillator  Strengths 

In  summary,  it  is  illustrative  to  list  the  oscillator  strengths 

used  in  the  present  calculations  compared  with  the  values  used  by 

,  ,  28,33 

other  authors. 


a  x  xxu 

REPORT 

28 

KECK,  ET  AL 

MEYEROTT 

NO  Beta 

.0015 

.006 

.  008 

NO  Gamma 

.0024 

.001 

.  0025 

N^  First  Positive 

.0095 

.025# 

.  02 

N^  Second  Positive 

.04 

.09 

.07 

N^+  First  Negative 

.04 

.  18 

.2 

O^  Schumann-Runge 

.  163 

.028 

.259 

33 


#A  revised  value  of  0.013  is  given  in  a  later  report  (Avco  Everett  Research 

Labs.,  Semi-Annual  Report,  June,  1961). 
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DISCUSSION  OF  RESULTS 


The  results  of  the  computations  are  presented  in  figure  form,  giving 
the  radiance,  either  spectral  or  total,  as  a  function  of  temperature  and 
density,  for  equilibrium  air  and  for  the  individual  radiating  species.  The 
data  is  given  in  abbreviated  form  in  Table  I,  which  includes  the  radiating 
system,  the  process  involved,  together  with  the  pertinent  critical  parameter, 
and  a  reference  to  the  techniques  used.  The  total  radiance  for  each  radiating 
system  is  given  by  an  equation  of  the  form 
J  =  A  (N)  exp  (-B/T) 

where  (N)  is  the  number  concentration  of  the  pertinent  species,  T  is  the 
.  o 

temperature  in  K,  and  A  and  B  are  constants.  The  equations  are  only 

valid  over  the  straight  line  portion  of  the  curves  in  the  cited  figures,  since 

they  ignore  self-absorption.  In  addition,  the  equations  for  the  ion  continua 

are  applicable  only  to  temperatures  above  8000°K. 

The  total  radiance  of  equilibrium  air  is  shown  in  Figures  1  through  10. 

The  first  three  of  these  figures  show  the  variation  of  the  radiance  as  a  function 

of  temperature  and  relative  density,  while  the  next  five  indicate  the  variation 

of  total  radiation  and  fraction  contributed  by  each  radiating  system  as  a 

function  of  temperature,  at  five  specific  relative  densities. 

Figure  9  compares  the  total  radiance  of  equilibrium  air  computed  by 

30 

our  techniques  with  the  earlier  results  of  Kivel  and  Bailey  .  The  spectral 
region  which  was  considered  by  Kivel  and  Bailey  is  between  .2  and  10M.  Our 
results  are  shown  for  the  region  .  2  to  10/Lt  and  also  for  the  extended  region 
.05  to  10/i.  It  is  obvious  that  this  additional  spectral  region  contributes 
significantly  to  the  radiation  from  air  and  cannot  be  neglected.  The 
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TABLE  1 


I 

! 


17 


♦Equations  are  valid  only  over  linear  region  of  curves  in  the  cited  figures 


comparison  of  the  data  between  .  2  and  10/i  shows  AVCO's  data  to  be  higher 

by  about  a  factor  of  2.  At  high  temperatures  the  difference  is  due  to  an 

30 

approximate  technique  used  by  Kivel  to  calculate  the  Kramers  radiation. 

,  •  29 

(Kivel  in  a  later  paper  indicates  that  the  effective  charge  assumed  was 
high  and  that  hence  the  results  are  high)  The  differences  at  lower  tem¬ 
peratures  are  due  to  different  theoretical  techniques  used  to  calculate 
radiation  from  molecular  bands  and  to  the  difference  in  the  electronic 
oscillator  strengths  used. 

Figure  75  predicts  the  radiant  intensity  of  equilibrium  air  using  normal 

53 

shock  temperatures  and  densities  as  given  by  Wittliff  .  These  predictions 
can  be  described  by  the  equation 

log  J  =  .  12  V  -  3.  5  x  10"ZH  +1.6 

where  V  -  ambient  velocity  (kFT.  /sec);  H  is  the  altitude  in  kFT.  and  J  is  the 

3 

radiant  intensity  in  watts/cm  -ster.  This  equation  is  accurate  for  velocities 
between  18  and  26  Kft/sec  and  altitudes  between  40  and  160  Kft.  At  higher 
altitudes  it  overestimates  the  radiation  by  a  factor  of  2. 

The  spectral  radiance  of  equilibrium  air  at  three  different  temperatures, 
and  at  a  relative  density  of  10  ^  is  given  in  Figures  11  and  12,  which  differ 
only  in  the  wavelength  range  covered. 

Figure  13  compares  the  spectral  radiance  of  air  obtained  in  our  com¬ 
putations  at  8000°K  and  a  relative  density  of  1  with  the  results  given  by  Keck 
28  3  3 

et  al  and  Meyerott  et  al  .  (Meyerott  gave  his  data  as  absorption  coef- 
ficients,  which  were  multiplied  by  the  blackbody  function  to  obtain  the  curve 
shown  in  Figure  13). 
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Figure  14  compares  our  spectral  absorption  coefficients  with  those  of 

Meyerott^^  at  12,  000°K,  the  highest  temperature  given  in  that  reference. 

The  agreement  is  generally  good,  the  smoother  shape  of  Meyerott's  curve 

being  due  to  the  fact  that  this  author  considered  fewer  absorption  edges  in 

his  computation  of  the  deionization  continua,  and  these  are  the  dominant 

radiators  under  these  conditions.  Quite  probably,  particularly  at  the  higher 

temperatures,  the  chief  uncertainty  in  the  radiation  due  to  equilibrium  air  is 

the  uncertainty  in  the  concentration  of  some  of  the  radiating  species. 

The  radiance  of  the  oxygen  Schumann-Runge  bands  is  given  in  Figures 

15  through  19,  and  these  figures  also  show  the  comparison  with  the  results 

of  Kivel  and  Bailey,  again  high  because  of  the  early  oscillator  strengths 

5 

values  they  used.  The  data  of  Bethke  used  in  these  computations  has  since 

43 

been  corroborated  by  the  work  of  Treanor  and  Wurster 

The  results  for  the  nitrogen  first  positive  system  are  given  in  Figures 

20  through  24.  Since  this  is  a  transition  between  two  excited  states,  the 

oscillator  strength  cannot  be  obtained  from  absorption  measurements  carried 

49 

out  at  low  temperature.  We  have  utilized  the  data  of  Wurster  to  compute 
an  "effective  oscillator  strength"  for  this  system  which,  matched  with  out 
computation  scheme,  reproduces  the  measured  absorption  coefficients. 

Similar  graphs  are  given  for  the  other  radiating  systems,  both  band 
systems  and  continua. 

52 

Recent  experimental  measurement  by  Wurster,  Treanor  &  Thompson 
have  indicated  an  additional  group  of  nitric  oxide  electronic  bands  located 
between  0.9  and  1.4  microns.  They  have  tentatively  been  identified  as 
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transitions  between  various  excited  state  of  NO,  and  Wurster  et  al  have 

indicated  that  they  can  be  roughly  approximated  as  having  an  equivalent 

excitation  energy  of  7  e.v.  and  an  effective  oscillator  strength  of  about  0.07. 

Measurements  in  air  at  6300°K  and  P/PQ  =  1  indicate  an  average  radiance  of 
3 

5  w/cm  -  steradian  -  micron  over  the  wavelength  interval  from  1  to  1.3 

microns.  The  absorption  coefficient  for  a  unit  density  of  nitric  oxide 

(p/p  )  =  1  is  about  0.  12  at  6300°K  for  the  previously  indicated  wave- 

o  NO 

length  interval  (1  -  1.3M). 

The  radiation  from  this  system  is  not  included  in  the  present  tabulation. 


20 


APPLICATION  OF  RESULTS 


The  results  of  the  computational  programs  are  presented  in  graphical 
form  to  permit  application  to  various  problems  involving  radiation  from 
air  or  its  component  gases. 

The  data  is  presented  insxch  a  form  that  it  can  be  conveniently  used 
to  determine  radiation  from  equilibrium  air,  from  the  individual  radiating 
species  in  equilibrium  air,  or  from  individual  radiating  species  in  non¬ 
equilibrium  air.  In  each  case,  one  need  only  find  the  pertinent  figure  (see 
Table  I)  and  read  off  the  radiationfor  the  appropriate  temperature  and  con¬ 
centration  (or  relative  density  -  as  the  case  may  be).  It  should  be  noted 
that  f-  r  non- equilibrium  compositions  the  pertinent  concentrations  are  those 
of  the  reactants,  and  caution  must  be  exercised  in  using  some  of  the  figures 
(particularly  Fig.  60  through  68)  which  for  computational  convenience  are 
given  in  terms  of  the  concentration  of  reaction  products.  Under  conditions 
of  chemical  equilibrium  the  two  concentrations  of  reactants  and  reaction 
products  are  related  by  the  equilibrium  constant,  and  hence  either  one  can 
be  used. 

The  quantity  tahulated  is 
, _  -KL.  _  _4 

(1-e  )  c  T  ,  2 

J_  =  - - -  watts /ster. -cm  (17) 

T  1 T 


where  L  =  1  cm 

^  T  „  <7T4  -KL.  ^  . 

For  J  «  -  ;  ( 1  -e  )  <  .  1 

t  It 


4 

KQ-  T 
It 


watts /ster.  -cm 


(18) 
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The  application  of  the  results  to  two  specific  problems  is  considered 
in  some  detail  in  the  next  two  sections: 

1.  Irradiance  from  a  Distant,  Optically  Thin  Source 

To  determine  the  radiation  emitted  by  a  known  volume  of  gas,  as 
viewed  at  a  distance,  one  simply  takes  the  appropriate  value  of  J  from 
the  tables,  multiplies  it  by  the  volume  (in  cubic  centimeters)  and  divides 
by  the  slant  range  (also  in  centimeters)  squared  to  directly  obtain 
irradiance  in  watts/centimeter^. 

This  computation  applies  equally  well  to  spectral  radiation  or  to 
total  radiation,  and  to  either  radiation  from  equilibrium  air,  or  to 
that  from  any  particular  radiating  system  of  air. 

The  only  assumption  made  is  that  the  radiating  layer  is  "optically 
thin".  This  is  generally  true  when  considering  gas  caps  associated 
with  re-entry  vehicles.  The  definition  for  an  optically  thin  radiator 
here  is  that  the  total  spectral  emissivity  K^L  is  less  than  about  0.  2 
or,  as  applied  to  air  radiation,  that  the  total  emissivity  K,j,L  is  less 
than  0.1.  (This  takes  into  account  that  in  high  temperature  air  the 
radiation  is  not  distributed  spectrally  in  the  same  manner  as  blackbody 
radiation).  In  general,  if  one  deals  with  radiating  layers  more  than  a 
few  centimeters  thick,  it  is  wise  to  check  if  the  gas  is  actually  optically 
thin  by  evaluating  KTL  (where  L  is  the  path  length  in  the  line  of  sight 
direction),  or  by  performing  the  equivalent  calculation;  i.  e. ,  taking 
the  product  of  (J^,)  (L)  and  dividing  by  J^,  for  a  blackbody  radiator  at 
the  same  temperature.  For  convenience,  curves  of  the  total  blackbody 
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radiation  (<7T  /iT)  or  a  convenient  fraction  thereof  are  included  in 
Figures  3  through  8.  A  typical  spectral  blackbody  distribution  curve 
for  8000°K  is  included  in  Figure  13. 

In  general,  the  plotted  data  is  designed  to  be  applied  only  to 
regions  of  constant  temperature  and  -  at  least  approximately  -  constant 
density.  The  hazards  of  using  an  "average"  temperature  are  best 
indicated  by  noting  the  dependence  of  the  total  air  radiation  on  temperature 
in  Figure  3. 

One  important  special  case  of  varying  temperature  and  density 

distribution  is  the  equilibrium  flow  field  around  a  hemispherical  body  in 

supersonic  flow.  Computations  were  performed  to  obtain  the  radiation 

from  the  bow  shock  region  of  hemispherical  bodies  (at  zero  angle  of 

attack)  at  several  different  altitudes  and  velocities.  It  was  found  possible 

to  relate  the  total  radiation  to  the  radiation  from  the  stagnation  region 

and  the  cube  of  the  radius  of  the  hemisphere  such  that 

3  , 

J.  =  constant  J.  R  watts/ steradian  (19) 

ii 
total  stagn 

3 

J  stagnation  (watts/cm  -  steradian)  is  a  function  of  T  stagnation  and 

p  .  ,  which  in  turn  are  a  function  only  of  velocity  and  altitude, 

stagnation 

"i"  refers  to  a  species  of  equilibrium  air,  or  to  the  total  equilibrium 
air.  The  following  constants  have  been  evaluated  for  hemispherical 
bodies: 
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Constants 


0.  25 
0.  36 
0.  28 
0.  21 
0.  17 
0.  31 
0.  28 
0.  61 


Oxygen  Atom  Free  -  Bound  Continuum  0.  23 

Nitrogen  Atom  Free  -  Free  Continuum  0.  19 

Oxygen  Atom  Free  -  Free  Continuum  0.  23 

The  constants  are  applicable  in  general  to  equilibrium  air 
radiation  from  bow  shock  regions  of  vehicles  with  hemispherical  noses, 
travelling  at  velocities  between  15,  000  and  25,  000  ft/sec.  and  at  altitudes 
between  250,  000  and  70,000  ft.  The  data  is  also  applicable  to  "blunt" 
vehicles,  of  the  sphere-cone  type,  for  cone  half  angles  of  less  than  15°. 

2.  Radiation  Heat  Transfer 

The  radiation  incident  or  one  square  centimeter  of  surface  from  an 


infinitesimal  thickness  of  an  optically  thin  layer  is  given  by 


2ffL  e 


watts /cm 


where  the  quantity  in  brackets  is  the  radiance  of  equilibrium  air,  as 
given  in  figures  1  through  8,  and  L  is  the  thickness  of  the  radiating  layer. 
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The  factor  two  is  exactly  applicable  only  for  the  case  of  an  infinite 
plane,  and  for  the  case  where  KL  (more  exactly  L  is  much  less  than 
1.  As  KL  approaches  1,  the  factor  2  also  approaches  1,  and  the  limiting 
case,  for  an  opaque  gas  is  simply  the  blackbody  equation: 


=  tr 
R 

32 


cr  T 

tT  i 


4  ,  2 

=  (7T  watts /cm 


(21) 


Mayer  has  considered  the  radiative  heat  transfer  for  typical  blunt 
re-entry  vehicles  travelling  at  velocities  of  the  order  of  20,  000  ft.  / sec. 
and  has  evaluated  that 


Qr  =  1 . 8  77  L 


ffT 


.  L 

watts  /cm 


(22) 
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Figure  1.  Total  Radiance  of  Equilibrium  Air  (3000  K  -  9000  K) 
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Figure  4.  Radiance  of  Equilibrium  Air  vs.  Temperature  (d/0q  =  10  ) 
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Figure  5.  Radiance  of  Equilibrium  Air  vs.  Temperature  (p/p  -  10  ) 
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Figure  7.  Radiance  of  Equilibrium  Air  vs.  Temperature  (0 /PQ  =  10°) 
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Figure  8.  Radiance  of  Equilibrium  Air  vs.  Temperature  (p/p  =  10  ) 
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Figure  10.  Comparison  of  Predictions:  Radiance  of  Equilibrium  Air  (18,  000°K) 
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Figure  12.  Spectral  Radiance  of  Equilibrium  Air  (1000  -  13,  000A) 
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Figure  15.  Radiance  of  Schumann -Rung e  Bands  (3000  K  -  9000  K) 
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Figure  16.  Radiance  of  O  Schumann-Runge  Bands  (10,000°K  -  25,  000°K) 
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Figure  21.  Radiance  of  First  Positive  Bands  (10,  000°K  -  25,  000  K) 
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Figure  26.  Radiance  of  Second  Positive  Bands  (10,000°K  -  25,  000°K) 


log  ([n2]l  KPARTICLES/CM2 ) 

■  _ T»  10000  •  K 


15000  ° K 
20000°  K 

25000  °  K 


r\i 


LOG  J  (WATTS /CM2 -STER.) 


2 


-I8h 


6  8  10  12  14 

LOG  ( [Nj]  L)  PARTICLES/CM2 

i.d6lo4id2  i.o^o1  io2  io"2  io°  io2 

"  ~  1  *  1  *  m  — -f — ■ - -yZ- - 


16 


18 


(T*K) 

4000 

6000 

8000 


10' 1  io°io'  IO2 


4 


Figure  31.  Radiance  of  N  +  First  Negative  Bands  (10,000  K  -  25,000  K) 
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Figure  38.  Spectral  Radiance  of  NO  Beta  Bands 
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Figure  39.  Radiance  of  NO  Gamma  Bands  (3000  K 
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Figure  42.  Spectral  Radiance  of  NO  Gamma  Bands 
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Figure  49.  Radiance  of  O  Free-Bound  Continuum  (0  Density  10 
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Figure  50.  Spectral  Radiance  of  O  Free-Bound  Continuum 
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Figure  51.  Radiance  of  N  Free  -  Free  (4000°K  -  9000°K) 
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Figure  52.  Radiance  of  N  Free  -  Free  (10,  000°K  -  25,  000°K) 
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Figure  53.  Radiance  of  N  Free  -  Free  (Cn]  L  e]  density  product  10 
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Figure  54.  Spectral  Radiance  of  N  Free  -  Free  Continuum 
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Figure  61.  Radiance  of  N  Deionization  (N  Density  10  -  10  ) 
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Figure  64.  Radiance  of  0+  Deionization  {0  Density  it)  ^  -  10^) 
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Figure  66,  Radiance  of  (N  +  O  )  Free  -  Free 
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Figure  67. 
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Figure  70.  Radiance  of  O  Deionization  (.05-.  2 AO  -  (O  Density  10 


( 1 0  “ 3  T)°K 


500  600  700  800  900  1000  1100  1200  1300 


LOG  J  (  WATTS  /  CM  d  -  STER. 


Figure  72.  Radiance  of  N  Deionization  (.05-. 2M)  -  (8000  K  -  25000  K) 
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